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I use quaternion free probability calculus — an extension of free probability to non-Hermitian 
matrices (which is introduced in a succinct but self-contained way) — to derive in the large-size limit 
the mean densities of the eigenvalues and singular values of sums of independent unitary random 
matrices, weighted by complex numbers. In the case of CUE summands, I write them in terms of 
two "master equations," which I then solve and numerically test in four specific cases. I conjecture a 
finite-size extension of these results, exploiting the complementary error function. I prove a central 
limit theorem, and its first sub-leading correction, for independent identically-distributed zero-drift 
unitary random matrices. 
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I. INTRODUCTION 



Model 



1. Definition of the model 



The main objective of this paper is to begin investigat- 
ing the following non-Hermitian random matrix model, 



Here, 



W = SP. 



S = wiUi + W2U2 



(1) 



(2) 



is a sum of L > 2 independent unitary random matrices 
of dimensions iV x TV, weighted by some arbitrary com- 
plex numbers wi, I — l,2,...,i. Everywhere, except 
subsection IHI Bl the U;'s will belong to the simplest cir- 
cular unitary ensemble (CUE). Moreover, 



P = A1A2 



(3) 



is a product oi K > 1 complex random matrices, 
where A/c, k — 1,2, . . . , K , is rectangular of dimensions 
Nk X Nk+i (hence, P has dimensions A^i x Nk+1: and 
there must be A'^ = iVi; the same are the dimensions 
of W), and where all the real and imaginary parts of 
the matrix elements of the A^ 's are independent random 
numbers with the Gaussian distribution of zero mean, 
i.e., concisely, 



where the ak 's are real positive parameters (which set the 
respective variances to be cr|/(2(iVfc7Vfc+i)-^/^)). Finally, 
any matrix entry of any U; is statistically independent 
from any entry of any A^. 

In this article, I will consider only the part S i.e., 
choose K ~ 0. The more general model W ([1]) will be 
left for a separate paper. 



2. Thermodynamic limit 

The tools I apply — quaternion free probability, with 
its quaternion addition law ([55)1 [H, 01, *-e-, an extension 
of the standard free probability addition law (|32t Q 
into the non-Hermitian realm — will allow to handle the 
above model only in the "thermodynamic limit," 



N = Ni,N2, 



Rk 



. . . ,Nk+i 00, 
= finite. 



(5) 



where the K finite parameters are called the "rectan- 
gularity ratios." 

However, in subsection IHI C[ I will conjecture a finite- 
size modification of the results for the model S (with 
cue's as the summands), featuring a simple form-factor 
(j67l) . which performs very well when numerically tested. 
The same form-factor should work for W. 
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(4) 



3. Mean densities of the eigenvalues and singular values 

I will be interested in the two simplest statistical prop- 
erties of the above model: 
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The mean density of the eigenvalues ( "mean spec- 



2 



tral density") of W (in this paper, only of S), 



1 ^ 



(6) 



There must be Nk+i = N (i.e., i?i = 1) for W to 
be a square matrix. Also, the averaging is per- 
formed with respect to the probability measure of 
W, and the complex Dirac delta is used because 
the eigenvalues \i are generically complex. 

• The mean density of the singular values, defined 
as the (real and non-negative) eigenvalues of the 
Hermitian random matrix H = W^^W, 



P-Si{x) 



K+l 



Nk + 1 



M^)) 



(7) 



(According to another terminology, these would be 
the singular values squared.) In this case Ri can 
be arbitrary, and H has dimensions Nk+i x Nk+i- 
The real Dirac delta is exploited since the singular 
values fj,i are real. 



More precisely, if one considers a bi-partite system 
consisting of a principal system A of size Ni and an en- 
vironment B of size N2, one may form a pure state as a 
linear combination of the product basis, 



Ni N2 

j=i 



(8) 



Now, a mixed state on the principal system is obtained 
by taking the "partial trace" over the environment. 



Tr(XXt) 



(9) 



Then, different probability distributions of the pure 
states (i.e., of X) lead to different ensembles of quan- 
tum density matrices. 

The model S ^ thus appears in the following applica- 
tion of the above program: Consider a bi-partite system 
{A,B) of size N X N (two "quNits"), and a maximally 
entangled state (the "generalized Bell state") on it. 



I* 



1 



AB/ 



N 

■E' 



|«M ® \VB- 



(10) 



B. Motivation 

1. Study of non-Hermitian random matrices 

The model W ([Ij is interesting from the mathemati- 
cal point of view, since it is non-Hermitian — and such 
random matrices have beautiful mathematical structure, 
more involved than Hermitian ones, plus multiple phys- 
ical applications, ranging from finances and biology to 
quantum physics (for a review, consult e.o., [5|). In par- 
ticular, the problems of summing {e.g., jll, [l,la-{l3l ) and 
multiplying {e.g., [T^ - [3^ ) non-Hermitian random ma- 
trices have been drawing considerable attention — and 
the model W includes both these operations. 



2. Applications to quantum entanglement 

The model W ([T]) arises in the theoryof random quan- 
tum states (see the textbook and js^, 113 for review; 
I base this introduction on these latter works). Such ob- 
jects are used for instance to describe states of a quantum 
system affected by noise, i.e., complicated interactions 
with an environment which can be regarded as random. 
Also, if one looks for generic properties of a complicated 
quantum state, one may assume it random. A random 
quantum state is defined by specifying a probability mea- 
sure in the space of density matrices p, i.e., Hermitian, 
weakly positive-definite (i.e., with non-negative eigen- 
values) and normalized (i.e., Trp — 1) matrices. One 
way to do this is to take any rectangular matrix model X, 
and then p = XX^/Tr(XX^^) is a proper random quan- 
tum density matrix. 



(A "typical" random state is entangled because there are 
much more entangled states than separable ones; the lat- 
ter form a set of measure zero in the set of all states.) 
Now, perform on this state L independent random local 
unitary transformations in the principal system A, 



\4>i)^{Vi®lN)\^\r,); 



(11) 



the resulting states remain maximally entangled. Form a 
probability mixture of these L states, i.e.^a superposition 
with coefhcients wi G [0, 1] such that J2i=i '^i — 1' 



\^k)=Y,wMi) = {S®\N)\'^' 



(12) 



1=1 



Finally, take the normalized partial trace over the envi- 
ronment B, which leads to the random mixed state © 
with X = S. 

The normalization, Tr(SSt) = A^(Ef=i HP) + ■ ■ •> 
where the dots are 2Y,i^v Re(u;;w;,Tr(U;U;,)) < N, if 
only the spectrum of the U; 's remains finite in the ther- 
modynamic limit (l5|). Hence, it is enough to study 
the model SS^, or practically equivalently (modulo zero 
modes), S^S, i.e., the singular values of S. 

To my knowledge, the mean density of the singular val- 
ues of S has been so far known only for all the wi 's equal, 
which is the "Kesten distribution" [s^; the correspond- 
ing mean spectral density is also known [l^, [ij] . 

The model P ([3]) appears in another setup: Consider a 
system consisting of an even number 2K of subsystems, 
with the following sizes, 



size A''! size N2 



■ , A2K-2,A2K-lj 

size Nii size Nk+i 



(13) 
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Consider an arbitrary product state 

I V'o) = \0)ai \0)a2 «) . . . ® \0)a2k ■ Now, form a 
random pure state by performing on ItJjq) independent 
random local unitary transformations acting on the 
following pairs of the subsystems, 

jlV-o). (14) 

By definition, the result is a product state with respect to 
this latter pairing, i.e., it can be expanded in the product 
basis as 

Afl N; Nk Nk + 1 

i^) = E E - E E 

n = l 12,4 = 1 !A',jif = 1 «ir + i = l {^^) 
■ \h)Ai «) 1^2)^2 ® |i2)-43 «) . . . \iK+l)A2K^ 

where the coefficients are collected into K matrices 
of rectangular dimensions Nk x A^fc+i. In the simplest 
case, they may be assumed Gaussian Consider fur- 
ther the maximally entangled states pUj) on the pairs 
(^^2,^3), • ■ • , iA2K-2,A2K-i), and perform the projec- 
tive measurement of {ip) onto the product of these Bell 
states, 

P ^ 1^, ® |*X.-2^2.-.)(*X.-2^2.-j) ® 

(16) 

which yields a random pure state describing the remain- 
ing two subsystems, Ai and A2K, 

!,/>)= 7'|V^)cx 

Ni Nk+1 ^^j^ 

E E i^]mK + i\h)AiE)\iK+l)A2K- 
il=l «K + l = l 

As a final step, take the normalized partial trace over 
A2K, thus obtaining the mixed state on Ai, ^ with 
X = P, whose statistical properties are directly related 
to the singular values of P. 

The mean density of the singular values of P, with the 
assumption of all the sizes Nk equal (i.e., all Rk = 1 (|5|)) . 
has been known as the "Fuss-Catalan distribution" [2J, 
HI]; the corresponding mean spectral density has been 
found in (29j . For arbitrary i?fc's, after a primary work 
on the K = 2 case [l^, both mean densities have been 
derived in [soj . 

The model W ([!]) arises when a combination of both 
the above procedures is applied, i.e., one takes a proba- 
bility mixture of L random pure states defined on the 2K 
subsystems, performs the projective measurement on the 
product basis of the {K — 1) maximally entangled states, 
and takes the partial trace over the last subsystem. 

According to my knowledge, the only known result con- 
cerns the mean density of the singular values of W for 
L = 2, with equal weights, and K = 1, which is the 



so-called "Bures distribution" [s^jHll. Therefore, gen- 
eralization to arbitrary L, K, weights, as well as consid- 
ering the eigenvalues — comprises an important research 
program, which will be accomplished in this and a forth- 
coming publication. 

3. Applications to random walks on regular trees 

I will just mention that the model S ©, in the CUE 
case, finds applications to random walks on L-regular 
trees [36j . 

C. Plan of the paper 

Section [TIJ "Quaternion free probability": I introduce 
the fundamental notions (Green function, M- 
transform, mean spectral density) of both Hermi- 
tian and non-Hermitian random matrix theory, in 
the latter case especially mentioning unitary matri- 
ces (jll Ap . I briefiy describe free probability theory 
of Voiculescu-Speicher and the addition algorithm, 
both in the Hermitian and non-Hermitian setting, 
the latter being the quaternion formalism, which is 
a crucial tool for this work (IIIBI) . I show how it 
is possible to reduce ( "rational Hermitization" ) the 
quaternion Green function for unitary matrices to 
the standard Green function (jll C|) . 

Section Hill "S umming free unitary random matrices 
— the eigenvalues" : I employ quaternion free prob- 
ability to the challenge of computing in the ther- 
modynamic limit the mean spectral density of a 
weighted sum of free CUE random matrices; be- 
sides a general ("master") equation, five specific 
examples are considered pil Ap . The central limit 
theorem for free identically-distributed zero-drift 
unitary random matrices, and its sub-leading term, 
are proven (IIII Bl) . I introduce and verify numeri- 
cally a modification of the mean spectral densities 
from subsection IIII AI which should be valid for fi- 
nite matrix dimensions and which uses the comple- 
mentary error function; I make a conjecture of its 
even broader application pil CI) . 

Section llVt "Summing free unitary random matrices — 
the singular values": I recall a conjecture relating 
the mean spectrum of the eigenvalues, when it is 
rotationally-symmetric around zero, and of the cor- 
responding singular values (jIV Ap . I exploit this hy- 
pothesis to find the master equation for the singular 
values, in the thermodynamic limit, of a weighted 
sum of free CUE random matrices; I also analyze 
the five examples from subsection IIII Al in this con- 
text (|rVB]) . 

Section |Vj "Conclusions": I summarize the results of 
the paper (jV Ap . and pose some unsolved problems 
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2. Unitary Green function 



complei; plane 



real asfis 



cuts of 
J, eigenvalues 



For unitary matrices U, the spectrum is not real, but is 
still one-dimensional (belongs to the centered unit circle 
C(0, 1), i.e., Xi = e'®% 0i G [0, 27r)), and the same formal- 
ism applies: The mean spectral density (in the variable 
9 G [0, 2n)) can be generically expanded in the Fourier se- 
ries, /3u(^) = 2¥(l + E„>i("i»e-'"«+m;e'"«)), where 
the coefficients m„ are called "moments." (The in- 
verse Fourier transform and triangle inequality imply 
|'7i„| < !•) Hence, the holomorphic Green function for 
U, 



FIG. 1; For Hermitian random matrices, the eigenvalues are 
real, and one needs a complex Green function of a complex 
variable. 



Gviz) 



2ir 



Pij{e)de 



1 + Mu(z) 



-M, 



u 



for \z\ > 1, 
, for \z\ < 1, 



(21) 

where the positive moments are gathered into a gener- 
ating function named the "holomorphic M-transform," 



dYH). 



Mu(z)^}_^^, for |z|>l. 

n>l 



(22) 



II. QUATERNION FREE PROBABILITY 
A. Hermitian and non— Hermitian Green functions 

1. Hermitian Green function 

"Quaternion free probability" is a version of 
Voiculescu's free probability calculus designed to handle 
non-Hermitian random matrices. To introduce it, let us 
however begin in the Hermitian realm. The "mean spec- 
tral density" of an x Hermitian random matrix H 
([7]) is conveniently encoded in terms of the "holomorphic 
Green function," 



Gu{z) 



1 

N 



Tr((zl.-H)-) = lf: 



1 



A,: 



(18) 

This is a meromorphic function, with poles coinciding 
with the mean (real) spectrum; in the large-iV limit, 
these poles coalesce into continuous intervals on the real 
axis, and this Green function turns into a holomorphic 
function on the whole complex plane except these cuts. 
The mean spectral density is retrieved from the holomor- 
phic Green function taken in the vicinity of these cuts, 
through 

Ph(A) = --^ lim (Gh(A + ie) - Gh(A - ie)) , (19) 

ZTTl e-i.O+ 

as follows from the representation of the real Dirac delta, 

S{X) = — 2^ liine_>o+ ( iTie ~ X^)' illustrated in 

figure [H Remark that alternatively to the Green func- 
tion, one often prefers the "holomorphic M-transform," 



3. Non-Hermitian Green function 

For non-Hermitian random matrices X, the ap- 
proach must be altogether different because the eigen- 
values are generically complex, and in the large-A'^ 
limit occupy on average some two-dimensional do- 
main T). The mean spectral density is now defined 
through the complex Dirac delta (|6]), whose represen- 
tation 5'^'^\\ — \i) = \d\' limg_j.o is at the roots 
of the definition of the "non-holomorphic Green func- 
tion" [liiil, 

Gx(z,z*) = hm lim — ( ) = 

^^o^^-^frt\|z-A.|' + e2/ 

= lim lim — Tr/- — — r- 

e^ON^ooN \(zlAr-X)(z*lAr-Xt)-he2liv 

(23) 

(with the convention for matrix division A/B = AB~^), 
since then the mean spectral density is obtained simply 
by taking a derivative, 

p^{z,z*)^-d,*G^{z,z*), for zeV. (24) 

TT 

(There are known intricacies concerning the order of lim- 
its in (1^51) . but I will not be bothered by them.) An 
equivalent object, often handier, is the "non-holomorphic 
M-transform," 



Mx(2,2*) = zG:i^(z,z*) - 1. 



Mh(z) = zGn{z) - 1. 



(20) 



(25) 



The non-holomorphic Green function (|23|) is a more 
complicated object than the holomorphic counterpart 
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B. Quaternion free probability in a nutshell 



quaternion space 




complex plane 



"Free probability," a theory initiated by Voiculescu 
and coworkers and Speicher Q , is a non-commutative 
probability theory (an instance of it being random ma- 
trix theory) endowed with a proper generalization of the 
classical notion of statistical independence, called "free- 
ness." Qualitatively, random matrices are free when not 
only are the entries of the distinct matrices statistically 
independent, but also when there is no angular corre- 
lation between them. I will not delve into details, just 
outline one important result of free probability, the "ad- 
dition algorithm." 



FIG. 2: For non-Hermitian random matrices, the eigenvalues 
are complex, and one needs a quaternion Green function of a 
quaternion variable. 

(|18p due to its denominator quadratic in X. Hence, a 
linearizing procedure has been proposed [lo| to introduce 
the "matrix-valued Green function," which is a 2 x 2 ma- 
trix function of a complex variable, 

gx{z,z*) = lim lim ^bTr/(Z,® ljv-XD)"'\, 

(26) 

where 

while bTr ("block-trace") turns a IN x 2A^ matrix into 
a 2 X 2 one by taking trace of its four N x N blocks, 

This need to leave the complex plane — into what 
we will see (paragraph IIIB 3p is the quaternion space 
— is shown in figure [2] Now, (|26)) is already linear 
in X, with the structure mimicking that of the holo- 
morphic Green function ([T51) . Its upper left element 
is precisely the non-holomorphic Green function (|23p. 

= Gx{z, z*). Its lower right element carries 
no new information, [Gxiz, z*)]22 = [Gx{z, z*)]l^. More- 
over, the negated product of the two off-diagonal ele- 
ments (being a non-negative real number), 

Cxiz.z*) = -[gxiz,z*)]i2[g^{z,z*)]2i, (29) 

has been shown [3, to describe correlations between 
left and right eigenvectors of X (a property I will not 
exploit), and also plays a role of an "order parameter": 
it is positive inside the mean spectral domain V and zero 
outside of V (because Cx{z, z*) cx e^, and the regulator 
e can be set to zero outside of V). From a practical point 
of view, once one finds Cx{z,z*) inside V, then setting 
it to zero and solving for z ^ x + iy yields an equation of 
the borderline dV in the Cartesian coordinates {x,y). 



1. Classical addition algorithm 

In classical probability theory, if two random num- 
bers Hi 2 are independent, then the PDF of their sum 
[Hi + H2) is derived using the "classical addition algo- 
rithm" : First, the PDF's of the constituents are encoded 
into the "characteristic functions," gHi_2{^) = (e'^^^'^), 
being complex functions of a real variable. Second, their 
logarithm is computed, rHi 2i^) = ^oggHi 2i^)- The in- 
dependence property then ensures that this object simply 
adds when summing the random variables, 

rHi+H2 {x) = rn^ (x) + r//^ (x). (30) 

Third, exponentiating the result leads to the character- 
istic function, carrying the full spectral information, of 
the sum. 



2. Hermitian addition algorithm 

Free probability provides a similar algorithm for Her- 
mitian random matrices; If Hi. 2 are free, then the mean 
spectral density of their sum (Hi -I- H2) is calculated as 
follows: First, one needs to have the holomorphic Green 
functions (|T8l) of the constituents. Second, one computes 
their "holomorphic Blue functions" (46)] , being functional 
inverses of the holomorphic Green functions, 

Gm,2 {Bu,,2 (z)) = 5hi.. {Gm,2 (z)) = z. (31) 

The freeness property then implies that these Blue func- 
tions obey 

Bii,+U2iz) = BuAz)+Bu2iz)-^. (32) 

Third, it remains to functionally invert the result to ob- 
tain the holomorphic Green function of the sum. (Let 
us mention that a more popular terminology is of the 
"i?-transform," RTi{z) = Bii{z) — 1/z, which is simply 
additive, i?Hi+H2(z) = Rhi{z) + Rh2{z)-) 



6 



3. Non-Hermitian (quaternion) addition algorithm 



"Rational Hermitization" procedure 



In [H, a straightforward extension to the non- 
Hermitian world has been proposed. The procedure hes 
on the observation that the matrix-valued Green func- 
tion looks analogously to the holomorphic Green 
function (jlSp considered in the vicinity of the mean eigen- 
value cuts, in the imaginary direction, i.e., Gh(A -I- ie) — 
while in order to compute the holomorphic Blue function 
PT|) . the holomorphic Green function must be known on 
the whole complex plane. Therefore, the authors of P, ^ 
felt compelled to replace the infinitesimally small regula- 
tor e in (j26p by an arbitrary complex number d (actually, 
real and non-negative d would be sufficient, since e is real 
and non-negative). 



1. Description of the method 

There exists a special class of random matrices which 
permit an explicit calculation of the quaternion Green 
function by means of a procedure called "rational 
Hermitization." To describe it, perform the matrix in- 
version in (IM)) . which gives the coefficients a, b of the 
quaternion Green function through the coefficients c, d 
of its argument. 



1 

'N' 



:Tr 



c*l 



N 



xt 



XXt - cXt - c*X + (|c|2 + |d|2) 1^, 



z le 
le z 



c id* 

id c* 



(33) 



(where c is also an arbitrary complex number), thereby 
defining the "quaternion Green function" as a quaternion 
function of a quaternion variable. 



a ib* 
ib a* 



. IbTr 

N 



(34) 

(I will henceforth refer to c, d and a, b as to the "coeffi- 
cients" of the respective quaternion.) 

The "quaternion addition algorithm" has then been 
proven: Let Xi 2 be free non-Hermitian random matri- 
ces. First, their quaternion Green functions (IMl) must be 
found (see subsection III CI for how it is done in the case 
of Hermitian or unitary matrices) . Second, these quater- 
nion Green functions are to be functionally inverted in 
the quaternion space, leading to the "quaternion Blue 
functions," 



(35) 



(38a) 



XXt - cXt - c*X + (|c|2 + |d|2) 1^ 



(38b) 



Consider a symmetry constraint on X such that X^ is 
a rational function of X. Two primary examples are 
Hermitian (H^ = H) and unitary (Ut = U""'^) matrices. 
Then on the RHS of P8ap - (j38bp one obtains rational 
functions of X, which in turn may be expanded into sim- 
ple fractions, and therefore written in terms of the holo- 
morphic Green function (jl8p of X. In other words, for 
this class of random matrices, knowing the holomorphic 
Green function is sufficient for knowing the (more com- 
plicated) quaternion Green function. 

The rational Hermitization procedure for Hermitian 
matrices has been outlined in [l|, [3| (it is needed e.g., 
when one adds a Hermitian to a non-Hermitian random 
matrix) . 



The freeness of the summands suffices to show that an 
analogue of ([5^ holds at the quaternion level, 



6xi+x.(Q) =Sxi(Q) + 6x.(Q)- Q"^ 



(36) 



Third, it is enough to functionally invert the result at the 
point -Ec=o to obtain the matrix-valued Green function 
(Uni) of the sum. 



X1+X2 



a ib* 
ib a* 



z 

z* 



(37) 



where a = Gxi+Xa (2, 2^*) and |5|2 = Cxi+Xs (2, 2*). 
(The regulator e may be set to zero here because it will 
be seen that the functional inversion to be performed 
in equation (1571) takes care by itself of regulating the 
singularities at the eigenvalues: (j37l) will always yield a 
"holomorphic solution," valid outside of V, and a "non- 
holomorphic solution," inside V.) 



2. Application to unitary random matrices 

For an arbitrary unitary random matrix w\J (where for 
further convenience I have included an arbitrary complex 
number w), one gets in this way the following coefficients 
of the quaternion Green function. 



1 

2^ 



\d\' 



& = --(l + Afu(w)+Mu(u)*): 
9 



i)Mu(u)*, (39a) 
(39b) 



where g = V((|c| - \w\)^ + |dp)((|c| + \w\)^ + |d|2) and 
u = f^d^wp + |cp -I- |d|2 -\- g), and recall the holomor- 
phic M-transform of U (f^. 
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III. SUMMING FREE UNITARY RANDOM 
MATRICES — THE EIGENVALUES 

After the self-contained introduction in sections |T] 
and ini I now come to the first main objective of this 
work, i.e., considering the model S ^ and exploiting 
the quaternion addition law p6|) . 



BsiQ) ^ B^.vAQ)' 



^B^,uAQ)-iL-m-\ (40) 



to calculate (j37p the non-holomorphic Green function 
of S, and consequently ([M]) its mean spectral den- 
sity. 



A. Summing free CUE random matrices 

Let us commence from solving this problem for all the 
U/'s belonging to the simplest unitary ensemble, the "cir- 
cular unitary ensemble" (CUE), defined to have all the 
moments zero, to„ = 0, i.e., Mu(u) — 0, i.e., the mean 
spectral density constant on C(0, 1), pcue(^) — l/2n. 



The master equation 



In this case, equations p9ap - (j39bp can be solved to 
provide the coefficients c, d of the quaternion Blue func- 
tion BwCUe{Q) through the coefficients a, b of the 
quaternion Q, 




Vl-4|wP|6p 



(41a) 



(41b) 



for some sign s = ±1. 

For the sum ([2]) of the CUE's with general weights wi, 
the quaternion addition law (j40|) along with (|37|) yield 
equations for the coefficients a, b of the matrix-valued 
Green function (^5]) of S, i.e., a = Gs{z, z*) (or better 
([25]), M EE Ms{z,z*) = za- 1) and \b\^ = Cs{z,z*), 



Z = Cl 



= di 



CL 



dL + 



(i- \)a* 
|a|2 + |6|2' 

{L-l)b 

H2+|6|2' 



(42a) 
(42b) 



where the q's and dis are given by (i4Tal) -(|4rb l) . 

This set (|ITa|) . (|41bp . (H^ . (|42bl) can be easily simpli- 
fied to the following form: M is a solution to the equation 



L + 2M = ^siJl + A\wi\ 



1=1 



M{M + 1) 



(43) 



for some proper choice of the signs si (yielding a mean- 
ingful solution), and then also 



Csiz,z*) 



M{M + 1) 



(44) 



This is the master equation for summing free CUE ran- 
dom matrices. 



Remark 1: The master equation p3|) and formula p4 
may be cast in a more elegant way. 



1=1 



where 

M{M + 1) 

Mi{Mi + 1) 



— C, / — 1, 2, . . . , L, 



(45) 

(46a) 
(46b) 



where for short C = Cs(2:,z*). This form is sur- 
prisingly symmetric: Equation (j46al) is identical as 
the set of definitions (ligb)) (c/. ([72l) - ((73bl) ). 



Remark 2: The above set (|45l) - (j46b[) also provides an 
algorithm for rewriting it as a single polynomial 
equation for M: (1) Find Ml from (gS]). (2) Sub- 
stitute it to (j46b|) . I — L. (3) Remove from it all 
the squares Mf, I = 1, . . . , L — 1, by using (j46bp . 
; = 1,...,L 1. (4) Find from it Ml-i (this 
requires solving a linear equation), and return to 
step (2), with I — L — 1. The procedure ends 
at Mq = M. One may also note that all the co- 
efficients of the resulting polynomial equation are 
themselves symmetric polynomials in the (L+l) ar- 
guments Z = 0, 1, . . . , L, where wq = z. How- 
ever, in all the five examples presented below, it is 
better not to use this algorithm, but rather succes- 
sively square the square roots in P5)) . as this leads 
to a polynomial equation of a lower order. For more 
general cases, however, the former algorithm seems 
inevitable. 



Remark 3: M depends only on 
R=\zl 



(47) 



i.e., the mean spectral density (I24p is rotationally- 
symmetric around zero, and may be expressed as 
ps{z,z*) = ^dji2M. Hence, we may focus on in- 
vestigating just its radial part. 



priR)^27:R ps{z,z*)\..j,^ 



dM 



(48) 



Remark 4: Only the absolute values, and not the com- 
plex arguments, of the weights wi are relevant. 
Hence, it is enough to study real and positive 
weights. 
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annulus of radii, 



FIG. 3: The 5 • lO'^ eigenvalues of the sum of L = 2 CUE ran- 
dom matrices of dimension A*' = 500, weighted with wi — 0.4 
and 102 = 1 — Wi = 0.6, obtained from 1, 000 Monte-Carlo it- 
erations (the points), and the theoretical annulus (|50a[) - (|50b|) 
(the dashed lines). 



Remark 5: As explained at the end of paragraph lll A 3[ 
once C is found, then C = is the equation of 
the borderline dV of the mean spectral domain. 
According to (|44|) (or (I46a|) ). this happens when 
M = or M = —1, i.e., in terms of the Green func- 
tion ((25| . G— 1/zorG^Q. This represents the 
matching on dV of the non-holomorphic quantities, 
valid inside T) {i.e., C, M, G), with the holomor- 
phic quantities, valid outside V {i.e., respectively, 
0, or —1, 1/z or 0). This means that the mean 
spectral domain D is either (1) a centered disk, en- 
closed by the circle R = i?ext — where the radius 
i?ext is found by substituting M = to the master 
equation — or (2) a centered annulus, confined by 
an external circle R — i?cxt {M = 0) and an inter- 
nal one R — Rint {M — —1). This is an instance 
of the "Feinberg-Zee single ring theorem" ana- 
lioj . which states that if the mean spectral density 
is rotationally-symmetric around zero, then V is 
either a disk or an annulus. 

I will now explicitly solve the master equation (I43p , or 
reduce it to a polynomial equation, in five cases. 



2. Example 1: L = 2 

As a first example of using the master equation 
consider a sum of L = 2 free CUE matrices, with arbi- 
trary weights wi^2- The master equation becomes linear, 
and yields the non-holomorphic solution, 



M 



1 



1 



i?2 



1 



(l"'l|-|"'2|)" 
i?2 



1 



(49) 



Both holomorphic solutions, M — and M — —1, are 
compatible with (|49p , which implies that I? is a centered 



Rcxt = Vklp + |W2|^ 



v\\ 



Wi\ 



\W2\ 



(50a) 
(50b) 



(It becomes a disk only when wi — W2, see figured (D).) 

The radial mean spectral density of S follows by ap- 
plying (|48l) to (|49t. 



pr{R) - 2R- 
I 



(51) 



1^21)^ 



1 



(h"ihh"2|)^ 



1 



for i?int < R < Rc-xt, and zero otherwise. (The apparent 
singularities lie outside of the annulus.) I show the exper- 
imental eigenvalues on the complex plane and mark the 
theoretical annulus in figure |3l Formula ()5ip is numeri- 
cally checked in figure H] (A) , finding perfect agreement 
in the bulk. 



3. Example 2: L = 3 

Consider now a sum oi L ~ 3 free CUE matrices, with 
arbitrary weights wi,2.3- The master equation can be 
transformed into a fourth-order polynomial equation, 



{R^ - vl) {R' - vl) {R^ - vl) {R^ - vl) 
+ Af3(^9i?s - 28Ai(i)i?^ + 10 (3/i(2) + 2/1(1,1)) R'^+ 
+ 12 (-/i(3) + /^(2,i) - 10M(i,i,i)) + vlvlvlvfj-V 
+ 2M^R^[lbR^ - 35/i(i)i?'' + {2b^i(2) + 14Ai(i,i)) R^ 

- 5/i(3) + 5/1(2,1) - 42/1(1,1,1)) 
+ AMR^ (lli?4 - 18/t(i)i?^ + 7/i(2) + 2/1(1,1)) + 
+ 24i?'' -^(1)) = 0, 



(52) 

where the basis of the monomial symmetric polynomials. 



7rePcrm(p) 7eSubs(L,#p) s=l 



27r(s) 



(53) 



where p is a partition, Perm(p) the set of all its permuta- 
tions, Subs(L, #p) the set of all the #p-element subsets 
of {1, . . . also Vl = I]fc=i(l - 2(5fe;)|wfc|. 

Inserting Af = or M = —1 into ((52|) leads to the val- 
ues of the radii of the circles enclosing the mean spectral 
domain D, 



Rint = Vmax(-yi,-F2,-F3,0), 



(54a) 
(54b) 
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MC histograms (solid)/theory (dashed), 
eigenvalues, 

ttiter= 1,000, N=500, L=2 




(Wi, W2) 

(0.4, «0.92) 
(0.4, 0.6) 




'MC histograms (solid)/theory (dashed), 
eigenvalues, 

|ttitei- 1,000, N=500, Li=3, 1^=7 




(Wl, W2) 


(0.2, «0.35) 


(0.2, i»0.06) 




FIG. 4: The radial part of the mean spectral density, ps"*(ii) (|48p . for (left to right, top to bottom): 

A: The sum of L = 2 CUE random matrices weighted with wi = 0.4 and (1) W2 = \/l — ~ 0.92, (2) W 2 

3 CUE ran dom matrice s weighted with (1) wi = 0.4, W2 
2 2" 



= 1 - lui = 0.6. 
0.69 {V is a disk). 



B: The sum of L = 3 CUE ran dom matrice s weighted with (1) w\ = 0.4, 11)2 = 0.6, W3 — \J\ — w\ — w\ 

(2) Wl = 0.3, W2 = 0.5, = \J\ — w\ — w\ ~ 0.81 (2? is an annulus), (3) w\ — 0.25, W2 = 0.35, 103 = 1 — wi — «;2 = 0.4 {V is 
a disk), (4) wi = 0.15, W2 = 0.25, ws = 1 — wi ~ 1^2 = 0.6 (O is an annulus). 

C: The sum of L = 10 CUE ran dom matrices such that ii = 3 weights equal wi, and L2 — 7 weights equal ■W2, where wi = 0.2 
and (1) W2 = V(l - Liwl)/L2 « 0.35, (2) w;2 = (1 - Liw^)/L2 « 0.06. 

D: The sum of L = 2, 3, 5, 10 CUE random matrices weighted with the same value w = I/v'L. One observes that as L grows, 
the bulk plots approach the straight line R — >■ 2R, corresponding to the GinUE. 

The theoretical results (|5ip . (|52p . (|55[1 . (|6ip . respectively (dashed lines) checked against numerical Monte-Carlo simulations, 
with matrix dimension N = 500, 1,000 iterations, and 100-bin histograms (solid lines). 



where VJ = ^ 2(5fe/)|wfep. One may check that 

the three numbers, —Vi^2,3, are either all negative or only 
one is positive; accordingly, I? is a disk or an annulus. 

I put forward a hypothesis that formulae analogous to 
(|54ap - (|54b|) hold for arbitrary L. 

The radial mean spectral density is given by (|48)) and 
equation (j52[) . and is numerically verified in figure 2] (B), 
with excellent concord in the bulk. 



4- Example 3: Arbitrary L and two "degenerate" weights 

Let now L be arbitrary, but the weights assume only 
two values — let Li weights be equal to some wi, and 
the other L2 weights to some W2 {Li + L2 ^ L). The 
master equation can be transformed into a third-order 



polynomial equation, 

M3 {R-L^\wi\ - ^2^21) {R + L^\wi\ - L2\W2\) ■ 
■ (i?- Lil^il +L2IW2I) [R + Li\wi\ + L2\w2\) + 

+ 2M^[lR^ - + 1) {Ll\u^i? + ^2^2!^) + 

+ {Ll\w,\^ - Ll\w2?fy 

/ (55) 
+ Mi^{L^ + L^L2) i?"*- 

- L ((2 + L2) + (2 + Li) L^|w2n + 

+ {L\\w,\^-Ll\w2?f) + 

+ LL1L2 - - L2\W2\^) = 0. 

Remark that ([55]) reduces to a quadratic equation when 
one of the lengths ii.2 equals 1, and further to a linear 
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equation when both lengths are 1, which case has been 
discussed in paragraph IIII A 21 

The values M = or Af = -1 substituted into ((55)) 
lead to the conclusion that if Li^2 > 2, then P is a disk 
of radius 



-Roxt = y/Li\wi\'^ + L2\W2\ 



(56) 



However, when one of the lengths equals 1, say Li, then 
V may happen to be an annulus of internal radius 



Rint = Vmax(|wi|2 - L2|w2p,0). 



(57) 



These expressions agree with the general hypothesis 
stated at the end of paragraph IIII A 3l 

The radial mean spectral density stemming from (j48p 
and equation ([55)) is numerically tested in figure S) (C), 
showing perfect coincidence in the bulk. 



5. Example \: Arbitrary L and equal weights. A central 
Itmit theorem 

Let again L be arbitrary, and this time all the weights 
be equal, wi = w. The master equation turns into a linear 
one, which yields the non-holomorphic solution. 



M = 



R2 _ fli' 



(58) 



where i?cxt is given by (|60)) . Let me also print the value 
of 



C 



(1 - i) - R') 



L J 



(59) 



The mean spectral domain T) is always a disk with 
radius 

i?cxt = (60) 
The radial mean spectral density (|^5)) . ()58p . 



pT{R) 



1 



R2_ 

L 



(61) 



for R < i?cxtj and zero otherwise. It is checked numeri- 
cally in figure 131(D), with perfect agreement in the bulk. 
This formula has been first derived in [l2| . and then in- 
dependently re-derived in [l3| . 

It is interesting to consider the limit L — >■ oo. It is clear 
from the above results that in order to arrive at a finite 
distribution, the weight w must scale as u; ~ Tak- 
ing L large, one finds that the mean spectral density (|CT)) 
becomes constant, ps{z, z*) — l/(7ri?^xJ, inside a disk of 
radius ~ which is the "Ginibre unitary ensemble" 
(GinUE) [53-113. I* is visible in figure S| (D). This is 
a first instance of the central limit theorem for unitary 
random matrices, here proven for free CUE's with equal 
weights; see paragraph IIII A 61 and subsection IIIIBI for 
generalizations . 



Example 5: L oo and small weights. A central limit 
theorem 



A simple generalization of the central limit theorem 
discussed in the previous paragraph would be to take 
L — >■ oo in the sum of free CUE's with arbitrary weights, 
only assuming \wi \ ^ 1. This may include situations such 
as wi = w{i/L)/y/L, where w{x) is a function on [0, 1]. 

Here it is convenient to use the master equation in the 
form (|i5)) - (|46bp . Since \wi\ < 1, (I46bp implies that also 
Ml ^ 1, if only we require a finite end result. Thus (|46bp 
becomes M; « —C\wi\-^. Denoting 



1=1 



(62) 



and supposing it to be finite and non-zero, the last equa- 
tion substituted into (|45)) gives — C ~ Af/cr^, upon which 
(|46ap finally yields the GinUE Green function G « z* ja^, 
i.e., the constant mean spectral density l/(7rf7^) inside a 
disk of radius i?cxt = cr- 



Central limit theorem 



Let me resume the discussion from paragraph IIII A31 
and again investigate the sum S ^ with equal weights 
wi = \/\/L, in the limit L — >■ cx), however now with- 
out any restriction on the probability distribution of the 
summands U;, except that it is identical for all the ma- 
trices and that its first moment (drift) vanishes, mi = 0. 
I have already shown that if the U; 's belong to the CUE 
class, the resulting mean spectral density is constant in- 
side C(0, 1) (GinUE with variance 1). I will now prove 
that for an arbitrary ensemble of the U/'s, the mean spec- 
tral density of S tends with increasing L to a universal 
distribution (independent of the details of the unitary en- 
semble, but dependent only on its second moment m2), 
namely the constant density inside a certain ellipse ()65l) . 

To do that, I substitute into the system of equations 
([MaP - pQb)) and (|i^ -(l42b | . along with (gg), the large- 
L expansions 



1 1 

a = ao + -^0.1 + ya2 + 

V i ^ 

b^bo + -^bi + y62 + ■ 

V L ^ 



(63a) 
(63b) 



and solve them order by order. Recall that a = Gs{z, z*) 
and = Cs{z,z*). 

At the leading order, one discovers 



ao 



1 



|m2| 



(64a) 



|5o|' = l + ^ ^ 



(l-\m2\' 



(64b) 



11 



Hence, the borderline of the mean spectral domain V 
(given by |6oP =0) is the ellipse 

((1 — Rem2)^ + (ImTO2)^) — 4(Imr7!,2).x?/+ 

'2 (65) 

+ ((1 + Rems)^ + (Imm2)2) = (l - |m2|^) , 

which has the semi-axes (l±|m2|), the angle 
ip G [— fif) from the x-axis to the major axis of 
the ellipse is i^s = iArg(m2) (with the convention that 
the principal argument lies in [— 7r,7r)), while its area 
reads 7r(l — |m2p). Furthermore, the mean spectral den- 
sity (|Ma|) . ([M)) is constant, ps{z, z*) = l/{n{l - |TO2p)), 
within this ellipse. This is the central limit theorem for 
free identically-distributed zero-drift unitary random 
matrices. 

I have also derived the next-to-leading-order expres- 
sions. 




(66b) 

(Note, |6p = |6oP + ;^26i6* + ....) They are propor- 
tional to the third moment 7713. 

C. Finite— size effects 

1. The erfc form-factor 

Random matrices exhibit freeness only in the limit of 
large matrix dimensions and consequently, all the 
above results hold only for A'^ — ^ c», so in particular 
are incapable of capturing the finite-iV behavior of the 
mean spectral density close to the borderline dV ( cf. the 
strong dumping visible in figure S]). 

However, this behavior has been described for a num- 
ber of non-Hermitian random matrix models {N x 
iV) which display rotationally-symmetric mean spec- 
trum and for which I? is a disk R = i?oxt — by 
a simple form-factor, ierfc((7(i? — i?ext) V^), where 
erfc(a;) = dt exp{—P) is the complementary error 

function, while q depends on the particular model. These 
models are: 

• GinUE [il, [13 (see also @); q has been derived 
explicitly. 

• P ([3]) with K — 2 and arbitrary R2 '28'] ; q has also 
been computed. 

• P with arbitrary K and arbitrary i?fe's [13]; the 
above form-factor has been conjectured and veri- 
fied numerically, with q treated as a parameter to 



be fitted by comparison with experimental data, as 
its exact form is yet unknown. 

• The "time-lagged covariance estimator" 
c oc ADA^^ , where A is a rectangular Gaus- 
sian random matrix and Dab = 5a+i,b [lH; q 
has been again treated as an adjustable parameter 
and remains to be derived. 



2. The erfc form-factor for weighted sums of free CUE 
random matrices 

Prompted by this performance of the erfc form-factor, 
I propose the following conjecture for the finite- A^ form 
of the mean spectral density of any weighted sum S ([2]) of 
CUE random matrices: The radial part of the N —i' 00 
formula for the density, p^s'^{R) (O, should be multi- 
plied by the factor of 

fN,,,,R,,s, (R) = ^erfc (qbSb (R - i?b) y^) , (67) 

for each centered circle R = Rh constituting a connected 
part of the borderline dT> (there can be only either one 
or two such circles, i.e., T) is either a disk or an annulus, 
as we know from the single ring theorem; recall Remark 
5 in paragraph IHI A l"|) . where the sign Sh is -1-1 for the 
external borderline and —1 for the internal borderline, 
while % is at the moment an adjustable parameter (one 
for each circle) to be found by fitting to experimental 
data, and eventually to be calculated. 

Figure [5] shows numerical tests of this conjecture for 
the same cases as in figure H] the proposal performs very 
well on all the examples. 



3. The erfc conjecture 

Since the erfc form-factor proves to perfectly repro- 
duce the Monte-Carlo data for the four very different 
matrix models described in paragraphs lHI C Tl and llll C 2[ 
I put forth a conjecture that (|67|) is valid for any non- 
Hermitian random matrix model whose mean spectrum 
possesses rotational symmetry around zero. 

Besides proving this hypothesis, another challenge is 
to express in each case the parameter(s) % through the 
parameters of the given model. 

In particular, it should work for W, provided that the 
U;'s belong to the CUE. 



IV. SUMMING FREE UNITARY RANDOM 
MATRICES — THE SINGULAR VALUES 

In this section, I will continue an analysis of a weighted 
sum of free unitary random matrices S ([2]) in the large- A?^ 
limit — focusing on the singular values. 
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0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 



FIG. 5: The same numerical histograms as in figure[4l checked against the finite-A'^ conjecture p^'^(R)fN,q,Re^t,+i{R) (when 
V is a disk) or Ps"^(i?)/]v,tji,i?ext,+i(-R)/jv,(j2,fli„t,-i(^) (when V is an annulus) ((67)l . where the least-squares values of the 
parameter(s) are: 

A: (gi.ga) ~ (1) (2.61,2.22), (2) (2.98,2.33). 

B: (1) g« 1.82, (2) {qi,q2) ~ (2.01,1.56), (3) g « 3.03, (4) {quqi) « (3.55,3.07). 

C:q^ (1) 1.49, (2) 4.27. 

D: q^ (1) 2.04, (2) 1.73, (3) 1.57, (4) 1.47. 



A. Conjecture about rotationally— symmetric 
spectra 

As noted in paragraph IIII A 11 Remark 3, the mean 
spectral density of a weighted sum of free CUE's (but 
not necessarily more general unitary ensembles!) is 
rotationally-symmetric around zero. In recent pa- 
pers [13, HH, the following conjecture has been claimed, 
which relates the mean eigenvalues and singular values 
of any non-Hermitian random matrix X with the men- 
tioned symmetry property: 

Step 1: The assumed symmetry of X can be re- 
stated as the rotational symmetry around zero 
of the non-holomorphic M-transform (P5|). 
Mx(z,z) = 9Jlx(i?^)- This allows to define its 
functional inverse, 

anx(OTx(^)) = z, (68) 

called the "rotationally-symmetric non- 
holomorphic iV-transform." 



Step 2: The random matrix X^X is Hermitian, thus 
one can always compute its holomorphic M- 
transform (pil)) . M^-i^{z), and the "holomorphic 
A^-transform" being its functional inverse, 

Mxtx(A^xtx(^)) (69) 

Step 3: The conjecture states that the two above A^- 
transforms remain in the following relationship, 

^xtx(^) = ^91x(^). (70) 



B. Summing free CUE random matrices 

1. The master equation 

It is straightforward to apply the hypothesis (|70)) to the 
master equation ()43|) : thus the master equation for the 
holomorphic M-transform, or better, the Green function 
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G = Gsts(-z); for free CUE's, reads 

L 



L-2 + 2zG = ^s,^l + 4|u;,|^zG2. (71) 



1=1 



Remark 1: (I7T|) may be recast as a set of equations 
analogous to (l45]) - (j46bp . with just one equation dif- 
ferent, 



1=1 



(72) 



where 



-zG^ = C, 



(73a) 

2 =G, ; = 1,2,...,L. (73b) 



Wl 



Remark 2: Recall that, as for any Hermitian random 
matrix, the Green function must behave at complex 
infinity z — >■ oo as G ^ 1/z. This will help us choose 
the proper solution of the master equation. 

Remark 3: It is a known fact [HI, that for any Her- 
mitian random matrix H, the end-points x* of 
the support of its mean spectral density, puix), 
are the branch points of the Green function, i.e., 
dG-H{z)/dz\z=xt = oo. Applied to the master equa- 
tion ([7T|) . this condition reads 



Si \Wl\ 



(74) 



1=1 Jl + A\wi\'x,Gl 



Thus the end-points arc found by solving the set of 
equations: ^ (with (z,G) = (a;*,G*)) and ((Til) , 
for x^, > 0. 

I will now revisit the five examples from subsec- 
tion 1111^ calculating the mean singular values. One will 
observe that in each case, the polynomial equation for G 
will have an order greater by 1 from the corresponding 
equation for the mean eigenvalues. 



2. Example 1: L = 2 

Consider a sum of L = 2 free CUE matrices, with 
arbitrary weights wi,2- The master equation becomes 
quadratic, and its solution with the proper large-z 
asymptotics reads 



Gsts(^) 



1 



\Jz- {\wi\ - \w2\f \Jz- {\wi\ + \w2\f 



(75) 



The mean spectral density of S^S follows from 
applied to ([75l), 

Psts(2^) = — , ^ 

7r^((|^i|-|^2|)'-x) (^x-i\w,\ + \w2\)') 

(76) 

for {\wi\ — \w2\)'^ < X < {\wi\ + |w2|)^, and zero other- 
wise. This is numerically checked in figure [H] (A), 
with perfect concord. Note that changing the argu- 
ment X = {\wi\ — 1^2!)^ -I- 4|wi||?«2|a;', where < < 1 , 
the density ([75]) becomes p'g\g{x') = l/{^T^yx'{l — x')), 
which is the "arcsine distribution." 



3. Example 2: L ^ 3 

For a sum of i = 3 free CUE matrices, with arbitrary 
weights wi,2,3j the master equation turns into a fifth- 
order polynomial equation, 

G'z {z ~ vl) {z - vl) {z ~ vl) {z - vl) + 
-I- 4G*z(^z^ — 3^(i)Z^ + (3^(2) + 2/X(i_i)) z— 

- M(3) + Ai(2a) - 10Ai(l,l,l))+ (-yy-j 

+ 2G3 (2z3 - 5Ai(i)z2 + 4/i(2)Z + ^1^21^3) - 

- 2G^ {z^ + VQViV2V:i) + 

-KG(-5z + 2^(i)) - 2^0, 

where recall paragraph IIII A 31 for notation. 

The mean density of the singular values is obtained 
by solving numerically (j77p at z = A ± ie, for small e 
([T^. which is compared with Monte-Carlo simulations 
in figure[6] (B), discovering perfect agreement. 

4. Example 3: Arbitrary L and two "degenerate" weights 

For arbitrary L — Li + L2, and Li weights wi and L2 
weights W2 , the master equation can be transformed into 
a fourth-order polynomial equation, 

GV (^Z-iL,\wi\+L2\w2\f) (^Z-{L,\Wi\-L2\w2\) 

+ 2G^z\L - 2) (z - L^lioip _ Ll\w2\^) + 

+ G^z^ (L^ - 6L + 6 + L1L2) z+ 

+ Ll{-LL2 + 2L-2)\wi\'^ + 

+ Ll (-LLi + 2L-2)|it;2p 

+ Gz{L - 2) (-2L + 2 + L1L2) - 
-(L-l)(Li-l)(L2-l)=0. 

(78) 

The mean density of the singular values numerically 
computed from ((78| . (IT9|) and Monte-Carlo simulated is 
shown in figure |6] (C), finding excellent agreement be- 
tween the two. 
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MC histograms (solid)/theory (dashed), | 
singular values, 

ttiter= 1 ,000, N=500, L=2 J 
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MC histograms (solid)/theory (dashed), 
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FIG. 6: The mean density of the singular values, pgig{x), for the same values of L and the weights as in figure |4l 



5. Example 4'- Arbitrary L and equal weights. A central 
Itmit theorem 



6. Example 5: L — >■ oo and small weights. A central limit 
theorem 



For arbitrary L and all the weights equal to some w, the 
master equation becomes quadratic, yielding one solution 
with the proper asymptotics at infinity. 



Gsts(2) 



1 Vf-(i-i)^Lt 



p2 

""-ext 



(79) 



where i?cxt is given by (pU]) 

Therefore, the mean density of the singular values ([T9| . 



'1 - 1] m i 



P2 _ ^ 



(80) 



for X e [0,4(1 — ■^)^oxt]' ^iid zero otherwise. This is the 
known "Kesten distribution" [s^, and has been derived 
in [l^ . It is successfully verified numerically in figure [6] 
(D). 

Taking the limit L — )■ cx) with i?cxt assumed fi- 
nite {i.e., w ~ leads of course to the Marcenko- 
Pastur distribution [Hi], psts(a;) = wTT (t - jw^V^^, 



for X £ [0, 4i?p^4] , and zero otherwise 



For L — >■ OO, in the presence of arbitrary weights 
such that \wi \ ^ 1, an analogous procedure as in para- 
graph IIVB 61 may be applied to the master equation in 
the form ([7^ - (j73bp : it gives a quadratic equation for the 
Green function, whose solution with the correct asymp- 
totic behavior at infinite z is 



Gsts(^) = 



1 



(81) 



where cr is defined in (15^ . and assumed finite and non- 
zero. This is again the Marcenko-Pastur distribution, 
with variance 



V. CONCLUSIONS 
A. Summary 

1. Main results 

In this article, I analyzed in the thermodynamic limit 
([5]) the weighted sum S Q of independent unitary ran- 
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doni matrices. The main results: 

• The "master equations" (|43| and dTl]) , which yield 
the mean densities of the eigenvalues and singu- 
lar values of S, respectively, in case of all the U;'s 
belonging to the CUE, but for arbitrary L and 
weights. 

• These master equations are either solved or trans- 
formed into a polynomial form, and numerically 
verified, in four cases: (1) i = 2 and arbitrary 
weights, (2) L — i and arbitrary weights, (3) any L 
and two "degenerate" weights, (4) any L and equal 
weights. The resulting polynomial equations have 
orders: For the eigenvalues: (1) 1, (2) 4, (3) 3, (4) 
1. For the singular values: (1) 2, (2) 5, (3) 4, (4) 2. 

• Two central limit theorems: (1) For the U/'s be- 
ing independent CUE matrices, and for arbitrary 
but small weights. The limiting distribution is the 
GinUE with variance (15^ . (2) More generally, for 
the Ui's being independent identically-distributed 
unitary random matrices of any probability distri- 
bution with first moment (drift) zero, and for equal 
weights, wi = l/^/Ti. The limiting distribution is 
an elliptic modification of the GinUE (|64ap - (j64b|) : 
I have also derived its first sub-leading correction 
(IMill - pbl) . 

• The conjecture about and numerical tests of the 
erfc form-factor ([67)) . 

These results, beyond being mathematically interest- 
ing, find applications e.g., in quantum entanglement the- 
ory and theory of random walks on regular trees. 

2. Method 

Another goal of this article was to advertise quater- 
nion free probability calculus (in particular, the quater- 
nion addition law as a conceptually simple, purely 
algebraic and efficient tool to compute mean spectral den- 
sities of sums of free non-Hermitian random matrices. 



B. Open problems 



two Hermitian random matrices. Hi = S^S and 
H2 = P^P. 

Step 2: The holomorphic M-transforms (|20p of both 
these matrices are known: Mhi {z) is given by equa- 
tion ([7T|l . while Mii^{z) has been found in [30| to 
obey a polynomial equation of order {K-\-l). Invert 
them functionally to obtain the respective holomor- 
phic iV-transforms ((69|) . 

Step 3: Use the "multiplication algorithm," well-known 
in free probability theory, valid for free Hi ^2, 

iVHiH.(2) = T^A^Hi(z)iVH,(z). 

Step 4: Invert functionally the result to get the holo- 
morphic M-transform M-<^^-<^{z), which contains 
the information about the mean density of the sin- 
gular values of W. 

Step 5: Since the mean spectral density of W must be 
rotationally-symmetric around zero (to be checked 
a posteriori), the conjecture (j70p finally yields the 
non-holomorphic M-transform M\\r{z,z*) (|25p . 
and consequently the mean spectral density of W 



2. Prove the conjectures 

One should also prove the mentioned conjectures: 

• The expressions for the external and internal radii, 
i?cxt and i?int, of the mean spectral domain V for 
an arbitrary weighted sum S of the CUE's (end of 
paragraph IIII A 3p . 

• The "erfc conjecture," concerning the valid- 
ity of the form-factor ((67)) . in the case of S 
with the CUE'S (paragraph IIIIC2p . and in full 
generality, i.e., for any non-Hermitian model 
with rotationally-symmetric mean spectrum (para- 
graph Unnil ■ 

• The conjecture relating the mean densities of the 
eigenvalues and singular values, in the thermody- 
namic limit, when the former displays rotational 
symmetry around zero (|70)) (subsection HVAl. 



1. Solve the model W 

In a forthcoming publication, I plan to repeat the con- 
siderations of sections IIIII (the mean spectral density) 
and IIVI (the mean singular values density) for the model 
W ([Ij , with the supposition that the U; 's belong to the 
CUE, and in the thermodynamic limit ([S]). The outline 
of the procedure: 

Step 1: Consider first the singular values of W, i.e., the 
model WtW = PtStSP. Through cycfic shifts, 
it is easy to relate this matrix to the product of 



3. Other open problems 

One might undertake the following research projects: 

• Investigate, using the quaternion addition law ()40p . 
the weighted sum S for the U; 's having more com- 
plicated probability distributions than CUE. In this 
case, generically, the mean spectral density will not 
be rotationally-symmetric around zero, hence, the 
conjecture (|70)) will not hold, and other means to 
reach the singular values will have to be devised. 
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• A major research program would be to search for 
the full JPDF of the eigenvalues of S (and of W) , 
for arbitrary L, K, weights, as well as — perhaps 
— other probability distributions of the Ui's and 
Afc's. 

Some other questions worth considering; 

• Calculate the "entanglement entropy" (see e.g.^ 
[M HI), 5x = -/dxpxtx(a;)^logx, for X = S 
and X = W. 

• Analyze more thoroughly the behavior of the mean 
singular values density Pxtx(2;), for X = S and 
X = W, near the end-points of the support (c/. 
Remark 3 in paragraph IIV B l]) . 
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